Search for: All records

Abstract We investigate the process of diffusive shock acceleration of particles with mass number to charge number ratiosA/Q > 1, e.g., partially ionized heavy ions. To this end, we introduce helium- and carbon-like ions at solar abundances into two-dimensional hybrid (kinetic ions–fluid electrons) simulations of nonrelativistic collisionless shocks. This study yields three main results: (1) Heavy ions are preferentially accelerated compared to hydrogen. For typical solar abundances, the energy transferred to accelerated helium ions is comparable to, or even exceeds, that of hydrogen, thereby enhancing the overall shock acceleration efficiency. (2) Accelerated helium ions contribute to magnetic field amplification, which increases the maximum attainable particle energy and steepens the spectra of accelerated particles. (3) The efficient acceleration of helium significantly enhances the production of hadronicγ-rays and neutrinos, likely dominating the one due to hydrogen. These effects should be taken into account, especially when modeling strong space and astrophysical shocks. 

Abstract A critical component of particle acceleration in astrophysical shocks is the nonresonant (Bell) instability, where the streaming of cosmic rays (CRs) leads to the amplification of magnetic fields necessary to scatter particles. In this work we use kinetic particle-in-cell simulations to investigate the high-CR-current regime, where the typical assumptions underlying the Bell instability break down. Despite being more strongly driven, significantly less magnetic field amplification is observed than in low-current cases, an effect due to the anisotropic heating that occurs in this regime. We also find that electron-scale modes, despite being the fastest growing, mostly lead to moderate electron heating and do not affect the late evolution or saturation of the instability. 

ABSTRACT The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $$\beta ^+$$ and $$\beta ^-$$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material. 

Abstract We present a survey of 1D kinetic particle-in-cell simulations of quasi-parallel nonrelativistic shocks to identify the environments favorable for electron acceleration. We explore an unprecedented range of shock speedsv_sh≈ 0.067–0.267c, Alfvén Mach numbers ${\mathcal{M}}_{\mathrm{A}}=5–40$, sonic Mach numbers ${\mathcal{M}}_{\mathrm{s}}=5–160$, as well as the proton-to-electron mass ratiosm_i/m_e= 16–1836. We find that high Alfvén Mach number shocks can channel a large fraction of their kinetic energy into nonthermal particles, self-sustaining magnetic turbulence and acceleration to larger and larger energies. The fraction of injected particles is ≲0.5% for electrons and ≈1% for protons, and the corresponding energy efficiencies are ≲2% and ≈10%, respectively. The extent of the nonthermal tail is sensitive to the Alfvén Mach number; when ${\mathcal{M}}_{\mathrm{A}}\lesssim 10$, the nonthermal electron distribution exhibits minimal growth beyond the average momentum of the downstream thermal protons, independently of the proton-to-electron mass ratio. Acceleration is slow for shocks with low sonic Mach numbers, yet nonthermal electrons still achieve momenta exceeding the downstream thermal proton momentum when the shock Alfvén Mach number is large enough. We provide simulation-based parameterizations of the transition from thermal to nonthermal distribution in the downstream (found at a momentum around $p_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}{v}_{\mathrm{sh}}\approx 3\sqrt{m_{\mathrm{i},\mathrm{e}}/m_{\mathrm{i}}}$), as well as the ratio of nonthermal electron to proton number density. The results are applicable to many different environments and are important for modeling shock-powered nonthermal radiation. 

Abstract Collisionless shocks tend to send charged particles into the upstream, driving electric currents through the plasma. Using kinetic particle-in-cell simulations, we investigate how the background thermal plasma neutralizes such currents in the upstream of quasi-parallel non-relativistic electron–proton shocks. We observe distinct processes in different regions: the far upstream, the shock precursor, and the shock foot. In the far upstream, the current is carried by nonthermal protons, which drive electrostatic modes and produce suprathermal electrons that move toward upstream infinity. Closer to the shock (in the precursor), both the current density and the momentum flux of the beam increase, which leads to electromagnetic streaming instabilities that contribute to the thermalization of suprathermal electrons. At the shock foot, these electrons are exposed to shock-reflected protons, resulting in a two-stream type instability. We analyze these processes and the resulting heating through particle tracking and controlled simulations. In particular, we show that the instability at the shock foot can make the effective thermal speed of electrons comparable to the drift speed of the reflected protons. These findings are important for understanding both the magnetic field amplification and the processes that may lead to the injection of suprathermal electrons into diffusive shock acceleration. 

Abstract The end of supernova remnant (SNR) evolution is characterized by a so-called “radiative” stage, in which efficient cooling of the hot bubble inside the forward shock slows expansion, leading to eventual shock breakup. Understanding SNR evolution at this stage is vital for predicting feedback in galaxies, since SNRs are expected to deposit their energy and momentum into the interstellar medium at the ends of their lives. A key prediction of SNR evolutionary models is the formation at the onset of the radiative stage of a cold, dense shell behind the forward shock. However, searches for these shells via their neutral hydrogen emission have had limited success. We instead introduce an independent observational signal of shell formation arising from the interaction between nonthermal particles accelerated by the SNR forward shock (cosmic rays) and the dense shell. Using a semi-analytic model of particle acceleration based on state-of-the-art simulations coupled with a high-resolution hydrodynamic model of SNR evolution, we predict the nonthermal emission that arises from this interaction. We demonstrate that the onset of the radiative stage leads to nonthermal signatures from radio to gamma rays, including radio and gamma-ray brightening by nearly 2 orders of magnitude. Such a signature may be detectable with current instruments, and will be resolvable with the next generation of gamma-ray telescopes (namely, the Cherenkov Telescope Array). 

Abstract The nonresonant streaming instability (Bell instability) plays a pivotal role in the acceleration and confinement of cosmic rays (CRs), yet the exact mechanism responsible for its saturation and the magnitude of the final amplified magnetic field have not been assessed from first principles. Using a survey of hybrid simulations (with kinetic ions and fluid electrons), we study the evolution of the Bell instability as a function of the parameters of the CR population. We find that at saturation, the magnetic pressure in the amplified field is comparable with the initial CR anisotropic pressure, rather than with the CR energy flux, as previously argued. These results provide a predictive prescription for the total magnetic field amplification expected in the many astrophysical environments where the Bell instability is important.